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Abstract 

Background, aim, and scope Matrix-based life cycle 
assessment (LCA) is part of the standard ingredients of 
modem LCA tools. An important aspect of matrix-based 
LCA that is straightforward to carry out, but that requires 
a careful mathematical handling, is the inclusion of 
sensitivity coefficients based on differentiating the 
matrix-based formulas. 

Materials and methods We briefly review the basic 
equations for LCA and the basic theory of sensitivity 
coefficients. 

Results We present the complete set of sensitivity coef¬ 
ficients from inventory to weighting through characteriza¬ 
tion and normalization. We show the specific formulas for 
perturbation analysis, uncertainty analysis, and key issue 
analysis. We also provide an example using the ecoinvent 
data. 

Discussion The limitations of the present approach include 
the restriction to small changes and uncertainties and the 
ignorance of correlation between input uncertainties. In 
contrast to common thinking, there is no restriction to 
normally distributed uncertainties: Every uncertainty distri¬ 
bution for which a variance can be defined can be 
submitted to the analytical uncertainty analysis. 
Conclusions This paper provides a useful set of tables for a 
number of purposes related to uncertainty and sensitivity 
analysis. 
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Recommendations and perspectives Although the formulas 
derived are not simple, they are straightforward to implement 
in software for LCA. Once this is done, the use of these 
formulas can become routine practice, enabling a key issue 
analysis and speeding up perturbation and uncertainty 
analysis. 

Keywords Derivatives • Life cycle interpretation • 
Matrix-based LCA • Sensitivity • Taylor series expansion • 
Uncertainty 

1 Introduction 

2 Background, aim, and scope 

By now, the matrix approach to life cycle assessment 
(LCA) has received wide recognition. Proposed in the early 
1990s (Moller 1992; Heijungs et al. 1992), it took a decade 
or so before this idea was embraced as a general accepted 
way of doing LCA (see Heijungs and Suh 2006 for a short 
history of the matrix approach to LCA). 

The matrix method has been developed to solve the 
inventory problem in LCA. With the inventory problem, we 
refer to the task of scaling all unit processes in the system in 
such a way that they exactly produce the reference flow (or 
functional unit) and to use this scaling such that the 
inventory table can be computed. Formulas such as 

g = BA"'f (1) 

(see Tables 1 and 2 for an explanation of the symbols 
involved) are now common in the specialized literature on 
life cycle inventory (LCI; Heijungs and Suh 2002; Peters 
2007; Tan 2008). 
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Table 1 Overview of symbols 
representing input data in the 

Symbol 

Name 

Dimension (rows x columns) 

Defined in 

LCA 

f 

Final demand vector 

Economic flows x 1 

Goal and scope definition 


A 

Technology matrix 

Economic flows x processes 

Inventory analysis 


B 

Intervention matrix 

Environmental flows x processes 

Inventory analysis 


Q 

Characterisation matrix 

Categories x environmental flows 

Impact assessment 


g 

Intervention totals a 

Environmental flows x 1 

Impact assessment 


h 

Category totals 3 

Categories x 1 

Impact assessment 

a See the note below on 

w 

Weighting factors 

1 x categories 

Impact assessment 

normalization 



The matrix approach has also been extended in various 
directions. For instance, research has been devoted to add 
details on: 

• the treatment of allocation and cutoff (Heijungs and 
Frischknecht 1998); 

• how to connect a process-based LCI to an input-output 
table (Suh and Huppes 2005); 

• how to efficiently compute an answer to the inventory 
problem (Peters 2007); 

• how to analyze the feedback structure of the system 
(Suh and Heijungs 2007); 

• how to calculate sensitivity coefficients (Heijungs 
1994). 

Most of these details indeed refer to exclusively inventory- 
oriented questions. The impact assessment follows the 
inventory results, so all issues related to cutoff, allocation, 
IO-based LCA, efficient algorithms, and structural path 
analysis are only interesting from an LCI point of view. For 
the last one mentioned, the sensitivity coefficients, this is 
different, however. 

Sensitivity coefficients are important for both uncertainty 
analysis and sensitivity analysis (Heijungs 1994). In the 
context of uncertainty analysis, they serve to establish 
essential information for a Taylor series expansion, and for 
sensitivity analysis, they provide the multipliers that enable 
one to distinguish sensitive from non-sensitive parameters, 
the so-called key issues for refined data collection (Heijungs 
1996). But uncertainty and sensitivity analyses are not only 
important in LCI but in impact assessment as well. 


Moreover, the impact assessment adds additional uncer¬ 
tainty to the already uncertain results of the LCI. Likewise, 
not only the sensitivity of inventory results is of interest 
but also (or perhaps even more so) the impact assessment 
results. 

The extension of the sensitivity coefficients from LCI to 
life cycle impact assessment was “left as an exercise” to the 
LCA practitioner. For instance, Heijungs and Suh (2002, 
p. 144) write that “In this way, all equations of LCA may be 
processed,” but they do not pursue this. It is a task that is 
not so often carried out, we guess, at least we have never 
seen the explicit results of such an exercise. This paper 
therefore aims to carry out this exercise and to make the 
results available. The formulas obtained can easily be 
implemented in matrix-based software for LCA. We have 
done so in CMLCA, a program for doing LCA, and some 
screenshots are shown at the end of the paper. 

In this paper, we first review the basic equations of LCA. 
Then, we proceed to review the theory of sensitivity 
coefficients in general and their form in LCI. This finally 
leads to a derivation and coherent presentation of the 
sensitivity coefficients for the entire LCA process. 

3 Materials and methods 

3.1 Basic equations for LCA 

The basic equations for LCA have been discussed in a 
consistent notation by Heijungs and Suh (2002). Below, we 


Table 2 Overview of symbols 
representing output results in the 

Symbol 

Name 

Dimension (rows x columns) 

Equation 

LCA and the formulas with 
which these results are obtained 

s 

Scaling factors 

processes x 1 

7 

< 

II 


g 

Inventory results 

Environmental flows x 1 

g = Bs 


h 

Characterization results 

Categories x 1 

h = Qg 


A 

Intensity matrix 

Environmental flows x 

A = BA 1 


h 

Normalization results 

economic flows 

Categories x 1 

W7=r 

hk 


h 

Category totals 3 

Categories x 1 

h = Qg 

a See the note below on 

W 

Weighted index 

lxl 

W = wh 

normalization 
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present two tables of symbols and equations connecting the 
fundamental concepts of LCA. 

With respect to normalization, there are two situations 
which require separate treatment. 

• A vector of intervention totals, g, can be defined for 
the reference situation. For instance, one can collect 
data on the emissions of C0 2 , S0 2 , etc., which then 
represent g 1? g 2 , etc. These then can be processed by 
the same characterization model to yield the vector of 
category totals, h. This then forms the basis of the 
normalization. Changing g will induce a change in h, 
but h itself will not be changed directly by the LCA 
practitioner. We will refer to this as normalization 
case 1. 

• Alternatively, the vector of category totals h can be 
known without a detailed specification of the underly¬ 
ing interventions. In that case, h is not an output result 
(in the sense of belonging to Table 2), but input data 
(in the sense of belonging to Table 1). Thus, h can be 
changed directly, and it will affect the normalization 
results and the weighted index. We will refer to this as 
normalization case 2. 


4 Results 


4.1 Sensitivity coefficients for LCA 


Table 2 provides a specification of the functions / in the 
context of LCA. So we are to insert the equations of 
Table 2 into the general equation for calculating a 
sensitivity coefficient (Eq. 4). Example calculations have 
been provided by Heijungs and Suh (2002, Eqs. (6.21), 
(6.26), and (6.29)). For the scaling factors, s, we have: 



and for the inventory results, g, we have: 



§t=< 6 » 

for the dependence on the elements of A, and 



Both approaches in fact appear in practice and are 
therefore elaborated below as separate cases. 

3.2 General theory of sensitivity coefficients 

There are various situations in which the stability of the results 
in terms of sensitivity for perturbations of the input data is of 
interest. In general, we may explore this issue as follows. 
Suppose that an output variable z depends on two input 
variables x and y and that the dependence is expressed by a 
function f. 

z =f(x,y). (2) 

The cmcial elements in determining sensitivity is the change 
of the result (Az) caused by a marginal change in x (Ax) and 
by a marginal change in y (Ay). This is expressed using the 
partial derivatives: 

dz dz 

Az = —Ax + — Ay. (3) 

dx ay 

Coefficients such as || and are referred to as sensitivity 
coefficients in the present context. Their evaluation requires 
a specification of the function f. Thus, with / specified, the 
sensitivity coefficient of/ with respect to input parameter x is 
defined as: 

dz_ = df{x,y) 

dx dx ' 1 ’ 


for the dependence on the elements of B, where 

1 if i = k represents the Kronecker delta. 

0 otherwise 

It is clear that the higher one moves in the sequence 
inventory-characterization-normalization-weighting, the 
larger the number of sensitivity coefficients there will be. 
Scaling factors only depend on the technology matrix, but 
the inventory results depend on the technology matrix and 
on the intervention matrix. 

Table 3 gives a tabular overview of the derivatives for all 
expressions of matrix-based LCA. 

One remark on these coefficients. Each of the formulas 
in Table 3 contains A” 1 or a symbol that depends on A -1 , 
such as s, A or h. Thus, we need to go through the process 
of matrix inversion to evaluate the sensitivity coefficients. 
Moreover, although Ciroth et al. (2004) argue that one can 
solve LCAs without calculating a full matrix inverse, 
Table 3 shows that we need the full matrix inverse when 
we wish to extend the analysis to include sensitivity and 
analytical uncertainty studies. 

4.2 Perturbation analysis 

Heijungs and Kleijn (2001) and Sakai and Yokoyama 
(2002) describe perturbation analysis as a way of investi¬ 
gating which input data are most decisive for the results in 
terms of their relative sensitivity. That is, given a system 
with the prototypical form 

^ =f(x,y ) (8) 
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Table 3 Overview of the sensitivity coefficients that express how LCA output results (columns) change if LCA input data (rows) data change 


ds k j 


Ogk/ 


dh k / 


dh k 


dW/ 


1 dai 


v 


! dh 


v 


'dq 


v 


(9g j (normalization case 1) 

'dhj (normalization case 2) 
'dwi 


-(A 

o 

o 

0 

0 

0 


ki s J 


^kiSj 


Sj&ik 


0 

0 

0 
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~ s j E ft/4 

i 

q ki^j 
gjdjk 

o 

0 

0 


~I k 


li 


tfkiSj 

hk 


gj_ _ h kgj 

hk ( h k y 

hqki 


8ik 


2 


ik 

hk 


0 


2 Oik 


-^/Ef Eft/4 

k k 1 
k V’t ) 

w (&-Mi 

V «' 

Wkh k qid 

t N 2 

Wihj 


hi 


07 

Example: row 3 and column 4 tell us that ^ = Cj^Sj 

Note that depending on the situation in normalization, either row 5 (• • • / dqy) or row 6 (• • • / dhj) applies 


one investigates dimensionless multipliers, such as 

dz x 
dx z 



The idea is that a small change of either of the input 
parameters (say in x, hence a change Ax) leads to a change 
in the result (z, hence Az), and that the relative change 
can be approximated by The results of Tables 2 and 3 
can be combined to yield a complete overview of these 
multipliers (see Table 4). 


4.3 Uncertainty analysis 


Although the Monte Carlo method is right now the most 
frequently applied method for studying the propagation of 
input uncertainties into output uncertainties (Lloyd and Ries 
2007), it has been recognized that this method may be too 
computationally intensive for application to large systems 
(Ciroth et al. 2004; Heijungs et al. 2005; Hong et al. 2008). 
Indeed, the present release of the ecoinvent data (v2.0, 
comprising almost 4,000 processes) has not been processed 
with Monte Carlo analysis, while the previous release 
(vl.3, 2,500 processes) contained such results, precisely for 
this reason. One way to address this is by means of smarter 
algorithms. Peters (2007) has proposed power series 
expansion, also in connection to Monte Carlo analysis. 
Moreover, there are more sophisticated sampling strategies 
than the naive Monte Carlo method, such as Latin 
hypercube sampling and response surface methods (e.g., 
Morgan and Henrion 1990). In this paper, we address it 
with a solution on the basis of analytical error propagation. 
The theory of analytical error propagation, using Taylor 
series expansion, has been proposed in LCA (Heijungs 
1994; Ciroth et al. 2004; Hong et al. 2008). Taylor series 
expansions are based on the approximation formula for 
calculating the variance of a stochastic result using 


Table 4 Overview of the relative sensitivity coefficients (multipliers) 
that express how small changes in input data propagate into changes 
in output results 


Multiplier 


Definition Formula 


^k(^ij) 

May) 

7 k(by) 

May) 

Mby) 

Mlij) 

ft(% 

M b v 


dsk/s k 


VkWij 


q k (gj) (normalization case 1) 
q k (hj) (normalization case 2) 

Aujj) 

u(bjj) 

Aqd 

co(gi) (normalization case 1) 
co(hi) (normalization case 2) 

u(wi) 


daij/atj 

dgk/gk 

day/ay 

dgk/gk 

dbij/bij 

dhk/hk 


daij/aij 

dhk/hk 


dbij/bij 

dhk/hk 


dqij/qij 

dh k /h k 

Ok! ij / Uy 

dh k /h k 

dby/by 

dh / h k 
dqy/qy 

dh k /h k 

dgi/gi 

dh k /h k 

dhi/hi 

dW/W 


daij/aij 
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Example: <r*(%) tells how many percent s k will change if a t j changes 1% 
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stochastic data (Bevington and Robinson 1994; Morgan 
and Henrion 1990). For a system of the form 

z=f(x,y) (10) 

it assumes the form 

var ^ = (lx) var( ^ + (iy) var ^ ( 11 ) 

df df . . 

+ 2— x — cov(x,y) 
ox oy 

where var(x) is the variance of the variable x and cov(x,y) 
represents the covariance between the stochastic variables x 
and y. 

In most cases, no covariance data are available, or the 
covariance can be assumed to be negligible as the 
uncertainties are in many cases independent. In those cases, 
we set covariance to zero and obtain 

var(z) = Q var(x) + (JQ var (y) (12) 

Such equations are for a general function f(x,y\ and they 
require the evaluation of the derivatives and |£. In the 
present case, we have specified the function (Table 2), and 
we have found equations for the derivatives (Table 3). So 
all ingredients are available to complete the structure of the 
uncertainty analysis. 

Heijungs and Suh (2002, Eqs. (6.73) and (8.87)) 
elaborate this only for the inventory vector, and then even 
only partially, and with a typo. The complete and correct 
expression for this is 


interpretation of these indices, we see summations over all 
economic flows and twice over all processes. Take the case 
of ecoinvent v2.0 where the number of processes and 
economic flows is almost 4,000; this variance is the result 
of almost 12,000 terms. Each of these terms is positive, so 

'-y 'y 

we may really consider {sjk ki ) var (ay) and (Sjk ki ) var (by) as 
the contribution that one individual var^) or var (by) 
makes to the total var(g>). We therefore define a number 
of dimensionless coefficients £, namely: 


G{gk,aij) 


{sjl h ) var(a, y ) 
var (g k ) 


and 


C(gk,by) 


(sj§i k ) 2 yar(bjj) 

var (g k ) 




as the relative contributions by each var (ay) and var (by) to 
the total var(g>). Naturally, 


Y^£{gk,aij)+Y^£(gk,bkj) = 1 (16) 

ij j 

so these indeed represent relative contributions to the 
variance of the total. 

Table 6 provides an overview of the expressions for the 
different types of key issues. 

4.4.1 Example 


var (g k ) = 22 ( s Aki) var(a,y) + 2J {sj) var (b kj ). (13) 

ij j 

The remaining expressions for the uncertainty of the 
impact assessment results are elaborated in Table 5. 

Again, mind that there are two expressions for var [hj^j 
and two expressions for var(JT) for both normalization 
variants. 

4.4 Key issue analysis 

Key issue analysis has been defined by Heijungs (1996) as 
the decomposition of the uncertainty of a result in terms of 
the contribution of the uncertainties of the input data. 
Morgan and Henrion (1990) refer to it as uncertainty 
importance, while Saltelli et al. (2000) use the (perhaps 
confusing) term sensitivity analysis. All expressions for the 
variance in Table 4 are the result of the weighted 
aggregation of a large number of variances of input data. 
For instance, the expression for var(g>) comprises two 
summations, one over i,j and one over j. Adding the 


We have programmed these equations into CMLCA 
(http://www.cmlca.eu/). As an example, we loaded the 
ecoinvent vl.3 (2,630 unit processes) and calculated the 
system for a reference flow of 1-kWh Swiss electricity, 
low voltage, at grid. In the example below, we have 
restricted the analysis to the inventory analysis and to the 
emission of Carbon dioxide, fossil, emitted to air, low 
population density. 

We performed a perturbation analysis, concentrating 
on the perturbations of the technology matrix. Figure 1 
shows a screenshot of the results, hence of the coefficients 
7 k (< 2 ij) for k is C0 2 and for all / (economic flows) and all j 
(processes). A not-too modern computer was able to 
complete the calculations within a few minutes. 

We also performed an uncertainty analysis. Figure 2 
shows a screenshot of these results, so tabulating var(g>) 
(along with some other statistics) for k is C0 2 . Again, 
results are obtained within a few minutes. 

We finally performed a key issue analysis. Figure 3 
shows a screenshot of these results, so tabulating C(gk> a ij) 
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Table 5 Overview of the variance of output results as a function of the variance of the input data 
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Table 6 Overview of the contributions to the variance of output 
results by the individual variance of the input data 

Contribution Formula 


<(%«/•/) 




C(§fo by) 


C(h h ciij) 


C (fib by) 


(( A ~7/7 var ( a ~ 

var(^) 

(sjki) var 

var fc) 


(sjS ik )\ai(bij) 
var (g k ) 



var (h k ) 


2 

(sjqq) var(%) 

™(h) 


C(b h qij) 


4 

{(hkibij) 
C (hk, tfij'j 
t(hk, &) 


(normalization case 1) 


(normalization case 2) 


CW«<r) 




C(W,g f ) (normalization case 1) 

£(W,hi) (normalization case 2) 

CWO 


(&) val A) 

var(Ajt) 



(A)’var(X) 




var(fiT) 



var(tT) 


Example: C(Sb a y) represents the contribution to the variance of by 
the variance of ay 


obtain a more reliable, less uncertain result by trying to find 
more accurate data for these four coefficients. 


5 Discussion 

Three important restrictions must be bom in mind. 

First, the formulas are for the “normal” LCI. This means 
that the equations become more complicated once we 
incorporate other features and developments, such as 
allocation and IO-based LCA. Heijungs et al. (2006) show 
how the formula for ^ changes when the inventory is done 
using a hybrid method, combining process-based and IO- 
based LCA. This modification may be propagated to the 
impact assessment level as well. 

Second, some of the formulas (namely those in Tables 4, 
5, and 6) are based on a first-order Taylor series 
approximation. That means that they are correct for small 
changes, uncertainties, and perturbations, but not neces¬ 
sarily for larger ones. To fix this, one may either include 
second-order terms (or even further). For instance, for the 
perturbations, we may improve by going from 


. dz 
Az = — Ax 
ox 



into 


dz 1 d 2 z 2 1 d 3 z 3 



This requires the evaluation of not only % but also of the 
second derivative or even beyond. Table 3 might be 

OX n 2 

extended to include expressions for ~ Sk , and similar. 

d { a ij) 

Alternatively, the approach by Sherman and Morrison 
(1950) can be used to provide an exact answer to this 
question. 

Third and finally, the formulas for uncertainty (Tables 5 and 
6) are based on the ignorance of the covariance between input 
variables. As we noted above, the complete first-order 
expression includes an additional term for the covariance 
between input data: 


and ((gk,by) for k is C0 2 and for all i (economic flows) 
and all j (processes). Again, a few minutes suffice. 

Here, we see that only four coefficients make up 80% of 
the uncertainty of the C0 2 emission. In other words, we can 


var(z) = (JQ var(x) + Q var (y) (19) 

df df . . 

+ 2— x — cov(x,j) 
ox oy 
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Fig. 1 Screenshot of CMLCA 
showing the results of a pertur¬ 
bation analysis, the sensitivity of 
an output result for small 
changes of the input data. Only 
multipliers smaller than -0.3 
and larger than 0.3 have been 
shown 


Perturbation analysis - analytical with first order Taylor expansion 


Alternative | [A1 ] Output of (G 25] etectiiaty, low voltage, at giidfCH ) ^ | 

1 Staling’ factors 2 Inventory analysis ] $ Chafacteiisation | Nasalisation | 5; Weighting | 



jnj xj 


Alternative 



Calculation— 

Resafciiate 


Settings ... 


Analyze 

Graph 



It is possible to carry out the program of this paper 
including covariance, but the expressions become much 
more complicated. For instance, for var(g>), we obtain 


var (g k ) = 22 (sjhi) var(a, y ) + (sf) var (b kj ) 


IJ 


J 


T 2 ^ ^ Sj?LkiS m ?iklZOv{aij, aim) 

iJJ.m 


+ 2 x y jS jSi^y{b kj M) 
hi 

+ 2 x E s j X ki sicov(a ij ,b k i). (20) 

iJJ 


The last three terms requires us to specify cov(a^a/ w ), 
co y(a k j,b kl ), and cov(a i j,b k i) for all combinations of 
processes, economic flows, and environmental flows. 
Although some uncertainties will be definitely correlated 
(for instance, the fuel input of a combustion process and 
the C0 2 emission of the same process may have a 
correlation close to 1), most uncertainties will be uncor¬ 
related, or any information on such correlations is lacking. 
The infrastructure needed (lots of memory for storing a 
number of covariance matrices, much more complicated 
formulas, much more data collection and estimation) will 
probably not offset the relatively limited gain of having a 
slightly more accurate computation. In the end, there is 
always something perverse about knowing the uncertainty 
with certainty. 

<£) Springer 


A question that always arises in connection to analytical 
error propagation is whether it only works for normally 
distributed uncertainties. This is not the case. The theory of 
analytical error propagation through Taylor series expansions 
(Morgan and Henrion 1990, p. 183 ff.) nowhere contains the 
assumption of normally distributed distributions. The only 
assumption is that for a first-order approximation, the 
function should be sufficiently close to a linear function 
within the range of uncertainty. This point has been addressed 
above. A practical issue is of course that the formulas require 
a specification of the variance, whereas most distributions are 
specified without an explicit variance. For instance, a uniform 
distribution is often specified in terms of its width and a 
lognormal distribution in terms of the geometric standard 
deviation (or its square, as in ecoinvent). Heijungs and 
Frischknecht (2005) provide formulas to easily calculate a 
variance from the standard parameters of a normal, lognor¬ 
mal, uniform, and triangular distribution. As we can see in 
Table 5, the propagated variances are the sum of a large 
number of terms. Following the central limit theorem (e.g., 
Morgan and Henrion 1990), such a propagated uncertainty 
will tend to become normally distributed provided that the 
input uncertainties are independent. 


6 Conclusions 

This paper has carried out the “exercise” that was left over by 
Heijungs and Suh (2002) and related work on deriving the 
complete set of sensitivity coefficients for matrix-based LCA. 
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Fig. 2 Screenshot of CMLCA 
showing the results of an 
analytical uncertainty analysis 
consisting of a number of 
characteristic parameters for the 
uncertainty of an output result 



7 Recommendations and perspectives 

For some, the formulas in Tables 4, 5 and 6 are intimidating 
and may offer little insight. But they are straightforward to 
implement in computer code. Once implemented, doing an 
uncertainty analysis is as easy as doing a Monte Carlo 
analysis: click a button and wait for the results. Likewise, 
doing a key issue analysis is as easy as doing a classical 
contribution analysis. 


We hope that the availability of these equations will 
stimulate developers of software, commercial or not, to 
implement analytical, Taylor series-based approaches 
toward uncertainty and sensitivity analysis. In particular 
for the key issue analysis, no good Monte Carlo 
approach is available, and the analytical solution using 
Table 6 provides an extremely powerful way of reducing 
the uncertainties in LCA. For perturbation and uncertainty 
analysis, numerical approaches can be used as well, but 



Fig. 3 Screenshot of CMLCA showing the results of a key issue analysis, the decomposition of the uncertainty of an output result in terms of 
contributions by the uncertainties of the input data. Contributions below 1% are not shown 
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these are extremely time-consuming for large LCA 
systems. 

As noted above, Monte Carlo analyses for large LCA 
systems may be unfeasible. This paper develops sensitivity 
coefficients that serve to derive the formulas for analytical 
error propagation based on first (or higher)-order Taylor 
series approximation. This has other limitations, for 
instance, relating to a more restricted range of validity. 
Therefore, we welcome the simultaneous development of 
more sophisticated sampling methods such as Latin 
hypercube sampling and response surface methods, as well 
as the development of more efficient algorithms such as 
those based on a power series expansion. Yet, even when 
this would overtake the analytical error propagation, we 
still see a role for the sensitivity coefficients in perturbation 
analysis and in key issue analysis. 

Open Access This article is distributed under the terms of the 
Creative Commons Attribution Noncommercial License which per¬ 
mits any noncommercial use, distribution, and reproduction in any 
medium, provided the original author(s) and source are credited. 
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